Chiral oscillations in terms of the zitterbewegung effect 
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Abstract 

We seek the immediate description of chiral oscillations in terms of the trembling motion de- 
scribed by the velocity (Dirac) operator ex. By taking into account the complete set of Dirac 
equation solutions which results in a free propagating Dirac wave packet composed by positive 
and negative frequency components, we report about the well-established zitterbewegung results 
and indicate how chiral oscillations can be expressed in terms of the well know quantum oscil- 
lating variables. We conclude with the interpretation of chiral oscillations as position very rapid 
oscillation projections onto the longitudinally decomposed direction of the motion. 

PACS numbers: 03.65.Pm, 11.30.Rd 



*Electronic address: alexeb@ifi.unicsunp.br 



1 



When Dirac had introduced his fundamental equation for relativistic quantum mechanics 
of spin one-half particles [1, 2] it became soon clear that in spite of its overwhelming suc- 
cess in describing very accurately the energy levels of the hydrogen atoms, this theory was 
replete of apparent inconsistencies as the presence of negative frequency solutions, the Klein 
paradox [3] and the zitterbewegung effect [4]. To be more specific, in first quantization, the 
relativistic Dirac equation predicts the existence for the spin one-half particle, under cer- 
tain conditions, of an oscillating time dependence in the average of position variable. This 
phenomenon, which is known as zitterbewegung, was formerly noticed by Schroedinger [4] 
as a consequence of the non commutative relation between the position operator x and the 
Hamiltonian H. The existence of the zitterbewegung effect is no stranger than the existence 
of negative frequency solutions since, in fact, such a trembling motion is only manifest for 
wave functions with significant interference between positive and negative frequency solu- 
tions of the Dirac equation. In atomic physics, the electron exhibits this violent quantum 
fluctuation in the position and becomes sensitive to an effective potential which could explain 
the Darwin term in the hydrogen atom [5] [23] . In the recent literature, we have witnessed 
some renowned interest in identifying the physical meaning of variables which coexist with 
the interference between positive and negative frequency solutions of the Dirac equation 
[6-9]. In this context, by following an almost identical line of reasoning as that applied 
for identifying rapid oscillations of the position, it is possible to verify that the expectation 
value of the Dirac chiral operator 7 s also presents a peculiar oscillatory behavior. In fact, 
a formalism with Dirac wave packets [5, 10] leads to the study of chiral oscillations [11]. 
In the standard model flavor-changing interactions, neutrinos with positive chirality are de- 
coupled from the neutrino absorbing charged weak currents [11] so that chiral coupled to 
flavor oscillations could lead to some small modifications to the standard flavor conversion 
formula [12]. Due to this tenuous relation between zitterbewegung and chiral oscillations, the 
question we shall answer in this manuscript is related to the immediate description of chiral 
oscillations in terms of the zitterbewegung motion, i. e. we will demonstrate that chiral 
oscillations are coupled with the zitterbewegung motion in a way that they cannot exist in- 
dependently each order. Moreover, if we observe that very rapid oscillations of the position 
variable can be decomposed into transversal and longitudinal polarization vector directions, 
we can also interpret chiral oscillations as position very rapid oscillation projections onto 
the momentum direction. The response to the posed question is obtained by means of very 
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simple mathematical manipulations, nevertheless, we believe that, in spite of its simplicity, 
it is important in the large context of quantum oscillation phenomena. 

All the following deductions concern with a free propagating fermionic particle, thus let 
us consider the covariant free particle Dirac equation given by 



{i^dn - m) ^(x) = 0, 



(1) 



where x = (t, x) and we have used the usual representation found in the literature [5, 13, 14] 
with c = h = 1. The plane wave solutions are expressed in terms of ip{x) = i]) + {x) + if)_{x) 
with 



ip + (x) = e' ipx ^ u(p) for positive frequencies and 



J+ipx] 



v(p) for negative frequencies, 



(2) 



where p is the relativistic quadrimomentum, p = (E,p) with E 2 = m 2 + p 2 , and the free 
propagating mass-eigenstate spinors are given by [16] 
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with u{p) (or v(p)) defined by u(p) = u\p)^° (or v(p) = v^{p)^°). We could also report to 
the Dirac wave packet formalism [13] where the general procedure consists in writing the 
general (complete) Dirac equation (wave packet) solution of Eq. (2) as 

V>(t,x) = / y^{b s (p)u s ( P ) exp [-iEt] + d* 3 (p)v s (p) exp [+iEt]}exp [ip- x] (4) 

where u(p) (or v(p)) is defined by u{p) = u\x)~f° (or v(p) = v^(p}^°) and p = (E,—p). By 
fixing the initial condition over ip(0,x) as the Fourier transform of the weight function 



ip(p - Po ) w = ^2{b s ( P )u s (p) + d*(p)v s (p)} 



(5) 



we get 



f dp r . 

W x) = J JT^y^iP - Po) ex P [IP-X\W 



(6) 



where w is some fixed normalized spinor. The coefficients b s (p) and d*(p) can thus be cal- 
culated by using the orthogonality properties of Dirac spinors. These coefficients carry an 
important physical information. For any initial state ip(0,x) given by Eq. (6), the negative 
frequency solution coefficient d*(p) necessarily provides a non-null contribution to the time 
evolving wave packet. This obliges us to take the complete set of Dirac equation solutions 
to construct the wave packet. Only if we consider a momentum distribution given by a 
delta function (plane wave limit) and suppose an initial spinor w being a positive energy 
mass-eigenstate with momentum p, the contribution due to the negative frequency solutions 
d*(p) will become null. 

By writing the free propagating particle Dirac Hamiltonian as 

H = a ■ p + (3m, (7) 

where ct = J2'l =1 o: k k = YH=ilolkk and (5 = 7„, we can easily verify whether a given 
observable O is a constant of the motion by means of the Heisenberg equation 

The Eq. (8) aids us to verify that, for instance, the free propagating particle momentum is 
a conserved quantity, 

^(p) = i{[H,p}) = 0. (9) 
Otherwise, the particle velocity given by 

d -{x) =i{[H,x]) = (a) (10) 

comes out as a non-null value. At first sight this appears quite reasonable [13] when we 
calculate the positive frequency solution (ip + (x)) averaged value 



(a)+(i) = Jd 3 x?pl{x)ot?p + (x) 



dfy p 
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3 = 1,2 

which represents exactly the expectation value of the relativistic velocity 2. Note, however, 
that the eigenvalue of a k is ±1 and corresponds to ±c, but we know that, for massive 
particles, the classical velocity cannot be equal to ±c. Moreover, once a k and a t do not 
commute when k ^ I, the measurement of the ^-component of the velocity is incompatible 



with the measurement of the y-component of the velocity. This may appear unusual since 
we know that p x and p y commute. In spite of this peculiarities, there is no contradiction 
with the above results. The plane wave solutions (3) which are momentum and energy 
eigenfunctions are not eigenfunctions of a^. Anyway the velocity related with a(i) does not 
represent a constant of the motion since 

j t {a) =i([H,a(t)]) =2i«p)-(aH)) (12) 

which can be regarded as a differential equation for a(t). Keeping in mind that p and H 
are constants of the motion we can easily solve the Eq. (12) in order to obtain [13] 

(a)(t) = (pH-i) + ((a(0) -pH- 1 )^- 2 " 1 *) (13) 

By following the same procedure, we observe that the spin angular momentum of a free 
Dirac particle represented by the operator S = 7 5 o: is not a constant of the motion since 

|(S) =<<[#,£]) = -2«a)xp) (14) 

With aid of the above result, it can be easily demonstrated that the particle helicity h = 
given by the projection of the spin angular momentum onto the momentum direction also 
is a constant of the motion, 

j t (h) = i([H,h}) = -((axp).p) = 0. (15) 

At the same time, the chiral operator 7 5 , which very often has its meaning confused with 
the helicity operator h, is not a constant of the motion [1 1] since 

j t (r>) = i([H^>]) = 2im<7V>- (16) 

The effective value of Eq. (16) appears only when both positive and negative frequencies are 
taken into account to compose a Dirac wave packet, i. e. the non-null mean value of (7 7 5 ) is 
revealed by the interference between Dirac equation solutions with opposite sign frequencies. 
It is important to note that an eigenstate of helicity can be read as an eigenstate of chirality 
only in the ultra-relativistic limit (mass m — 0) [16]. At time t — 0, the coefficients b s (p) 
and d*(p) used in the construction of the Dirac wave packet tp(x) can be chosen to provide 
a negative (positive) chiralitly eigenstate [12, 15] or, in the same way, to provide a helicity 
eigenstate (when hu s (p)(v 3 (p)) = ±^u s (p)(v s (p))). Once we have assumed that the initial 



chiral eigenstate[24] ip(0,x) is not only a superposition of momentum eigenstantes weighted 
by a momentum distribution of Eq. (5) centered around p , but also a helicity (constant of 
the motion) eigenstate obtained through the production of a spin-polarized particle, which 
formally occurs when we assume that the constant spinor w in the wave packet expression (6) 
is a simultaneous eigenspinor of 7 5 and h[25], we can make use of the following decomposition, 

j t (l*){t) = j t (cc-p(V.p))(t)= (L( a .p)(t)) (V.p)(t) + ( a .p)(t) (^(£.p>(i))(17) 

from which, when we substitute the Eqs. (9) and (16), we can establish a subtle relation 
between the chirality operator 7B (t) and the velocity operator a(t) in the following way 

!<7.X*) = (2h) (| <«)(*))• p. (18) 

The time evolution of 7 5 presents an oscillating character which can be interpreted as a direct 
and natural consequence of the oscillating character of the trembling motion described by 
a(t). The coupling of Eqs. (12) and (18) leads to the explicit dependence of 7B (t) on the 
momentum variable, ac(t), 

^( 7B )(t) = 4i/i (p-{a(t)H))-p. (19) 

Since (h), p and H are constants of the motion, we conclude that there will not be chiral 
oscillations without the "quivering motion" of the position. The constraint between the 
chirality operator 7 5 (i) and the velocity operator a(t) becomes more interesting when we 
take into account the complete expression for the current density ip(x)^ /Ji ip(x) (which leads 
to the averaged value of a(t)) and apply the Gordon decomposition [16], 

^(xM(x) = -^[(d^(x))iP(x)-^(x)(d^(x))] 

+ 2^0" $(x)a^(x)) , (20) 

where a^ v = | [7^,71/]. The spatial integration of the vector components of Eq. (20) gives us 

J tfic^a^ = ^ Jd A x{-i p(V^) - (V^) V] 

+ [V x pS^) - id t i^ctif)] } (21) 

where we have suppressed the x dependence. By using the Dirac equation solution expressed 
by (4), we can write the decomposed components of (a) as 

Jd^V x = 0, (22) 
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with p = (E, —p), where we have assumed a s = i] s > a ■ p r] 3 = (— l) s+1 5 ss ', and 

__L jtket&n,) = / j ^{E(3«.*[»:w:* l+I,E " - iw».w^ iH| ] 

+ 5>,['>>><e<^>e [+2 ' £ ' , - d.W6,( I i)e[- 2 " ! 'l]}(24) 

where n^j = 1 ± i2 when p = 3 and the unitary vectors 1, 2 and 3 correspond to mutually 
orthogonal directions. 

The Eq. (22) allows us to verify that the zitterbewegung does not get a contribution from 
the intrinsic spin dependent (XI) magnetic moment component which couples with external 
magnetic fields B(x). In fact, the zitterbewegung originates from the current strictly related 
to the internal electric moment. To clear up this point, it is convenient to consider the mod- 
ern and more precise interpretation of the zitterbewegung [17, 18]. From such a theoretical 
perspective, the zitterbewegung for a Dirac particle is clearly related to the separation be- 
tween the center of mass, which is related to the Foldy-Wouthuysen position operator, and 
the center of charge, which corresponds to Dirac's position operator x. It is the particle's 
charge at position x that is moving at the speed of light in circles of radius h/2mc around the 
center of mass, such that the average value of this velocity is related to the linear momentum 
of the particle. It is this separation, and thus the existence of an electric dipole moment 
with respect to the center of mass, which justifies the above statement. For both particles 
and antiparticles, then in pure positive or negative energy states, this internal motion of the 
charge around the center of mass exists but there are no chiral oscillations, which is no longer 
obvious. In the wave packet construction, it should be equivalent to assume only positive 
(or only negative) Dirac equation frequency solutions (like in ()11)) in the calculation of 
(observable) averaged values, in particular, for the 7 s time-derivative operator (which would 
lead to d(Y°)/dt = 0). However, the production of a fermionic particle as a chiral eigenstate 
(for instance, a neutrino), imposes the wave packet space-time evolution described by the 



7 



superposition of positive and negative frequency solutions of the Dirac equation [10]. It 
recovers the possibility of chiral oscillations for massive fermionic particles. In this context, 
the coupling between chiral oscillations and the zitterbwegung is recovered. 

By taking into account the (momentum direction components of the) Eqs. (23-24), we can 
turn back to the Eq. (19), which carries the main idea of this manuscript, and observe that 
chiral oscillations can be essentially constructed in terms of the longitudinal components of 
(a). By calculating the mean value of (ocH) and projecting it onto the momentum direction 
p, we obtain 

+ Ef «• [(m(p)d*Me^ - (E)d s ( P )bMe [ ^ Et] ]\ 

3 = 1,2 ) 

= \P\-jj^y m X> [&>K(iW +2iBtl - d.{p)bMel- 2iE *] , (25) 
which can be substituted in (19) in order to give 

j t (l,){t) = |^^£(2/iaJ(2itf)[^^ (26) 

in a way that the time evolution of the chirality operator could be written as 

(7.X*) = (7s)(0)+ /tSi| E( 2/ ^) [<*.(P)6.(P) (e [2tEt] -l)+h.c] 

= (7.X0) + j^y | £ [dM b.m (e^ - 1) + h.c] . (27) 

since 2ha s is equal to the unity for well defined spin up/down (helicity) states. 

The physical significance of the above results is discussed in a very subtle way in [11]. 
They demonstrated that, for a gaussian initial wave function and if the initial state has 
average chirality zero, the oscillations of the left-right (L-R) chiralities cancel and there is 
again no overall oscillation. It could be the origin of an apparent paradox. In order to 
reproduce the ideas presented in [11] we observe that for any mass-eigenstate represented 
by plane wave solutions of the Dirac equation with mass, the rest-frame wave function is 
always an equal mixture of both chiralities. This is easily verified when 

1 - t b 1 + r 5 

i> = + = *Pl + fpRi (28) 
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where ipL,R correspond respectively to chirality quantum numbers =)=1. It could be shown 
that, in the rest frame of a particle, we have 

W=|lM a . (29) 

Note that this result is not Lorentz invariant since a Lorentz boost is not a unitary transfor- 
mation. Thus, while the cross section is Lorentz invariant the chiral probabilities are not. 
This seems to suggest that probabilities measurements are chiral independent. We seem 
to have an argument against the physical significance of chiral oscillations. The reply to 
this objection based upon the Lorentz invariance is simply that in any given Lorentz frame 
chiral oscillations are manifestly important because of the chiral projection form (V-A) of 
the charged weak currents. The chiral probability variations produced by Lorentz transfor- 
mations (even if 7 5 commutes with the Lorentz generators) are automatically compensated 
by the wave function normalization conditions [26] and the Lorentz transformations of the 
intermediate vector bosons and other participating particles. 

In addition, the above discussion introduces the primary tools for accurately deriving 
the expression for the neutrino spin-flipping in magnetic field which can be related to chiral 
oscillations in the limit of a massless particle (ultra- relativistic limit). By correctly differing 
the concepts of helicity and chirality, we can determine the origin and the influence of 
chiral oscillations and spin-flipping in the complete flavor conversion formula. In some 
previous manuscripts [10, 12, 15] we have also confirmed that the fermionic character of 
the particles modify the standard oscillation probability which is previously obtained by 
implicitly assuming a scalar nature of the mass-eigenstates. Strictly speaking, we have 
obtained the term of very high oscillation frequency depending on the sum of energies in 
the new oscillation probability formula which, in case of Dirac wave-packets, represents 
modifications that introduce correction factors which, under the current experimental point 
of view, are not effective in the UR limit of propagating neutrinos in vacuum but which 
deserves, at least, a careful investigation for neutrino moving in the background matter where 
spin/chiral effects become more relevant [19, 20]. The physical consequences in environments 
such as supernova can be theoretically studied [21]. For instance, it was observed that 
neutrinos propagating in matter achieve an effective electromagnetic vertex which affects 
the flavor conversion process in a framework where preserving chirality can be established 
[22]. 
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Finally, just as a remark about this connection with neutrino physics, it is also to be 
noted that in this kind of analysis we have to assume that neutrinos are Dirac particles, 
thus making the positive-chiral component sterile. If neutrinos are Majorana particles [19], 
they cannot have a magnetic moment, obviating the spin- flipping via magnetic field interac- 
tions but still allowing the (vacuum) chiral conversion possibility via very rapid oscillations 
(zitterbewegung) . 

To conclude, we have shown how chiral oscillations can be mathematically related to 
the zitterbewegung motion described by the velocity (Dirac) operator ex. Once we have 
assumed the neutrino electroweak interactions at the source and detector are (left) chiral 
(■07 M (1 — ^ipWn) , only the component with negative chirality contributes to the propaga- 
tion, therefore, chiral oscillations can take place in the context of neutrino quantum oscil- 
lation phenomena. We remind that, in the usual (standard) treatment of vacuum neutrino 
oscillations, the use of scalar mass-eigenstate wave packets composed exclusively by posi- 
tive frequency plane-wave solutions is implicitly assumed. Although the standard oscillation 
formula can predict the correct result when properly interpreted, a satisfactory description 
of fermionic (spin one-half) particles requires the use of the Dirac equation as evolution 
equation for the propagating mass-eigenstates. Consequently, the the spinorial form and 
the interference between positive and negative frequency components of the mass-eigenstate 
wave packets leads to the possibility of chiral coupled with flavor oscillations which effectively 
introduce some small modifications to the standard flavor conversion formula [10] when it 
concerns with non-relativistic neutrinos. 
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